Persistent spin oscillations in a spin-orbit-coupled superconductor 
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Quasi-two-dimensional superconductors with tunable spin-orbit coupling are very interesting sys- 
tems with properties that are also potentially useful for applications. In this Letter we demonstrate 
that these systems exhibit undamped collective spin oscillations that can be excited by the appli- 
cation of a supercurrent. We propose to use these collective excitations to realize persistent spin 
oscillators operating in the frequency range of 10 GHz — 1 THz. 
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Introduction. — Spin-orbit-coupled two-dimensional 
(2D) electron gases (EGs) are the focus of great interest 
in the field of semiconductor spintronicJ^. This interest 
has been largely fueled by the hope to realize the vi- 
sionary Datta-Das "spin transistor"^ in which the on/off 
state is achieved by purely-electrical control of the elec- 
tron's spin in a spin-orbit-coupled semiconductor channel 
placed between ferromagnetic leads. Research in spin- 
orbit-coupled 2DEGs has been recently revitalized by 
theoretical^ and experimentaP studies of the spin Hall 
effect, in which a current traversing the sample generates 
a spin-current in the orthogonal direction. 

The study of the interplay between spin-orbit coupling 
(SOC) and superconductivity in 2D systems, stemming 
from the seminal works of Edelstein-^ and Gor'kov and 
Rashbsl^, has also gained impetu^l. There is a large vari- 
ety of systems in which SOC and superconductivity coex- 
ist: two examples of great current interest are i) 2DEGs 
in InAs or GaAs semiconductor heterostructures that are 
proximized by ordinary s-wave superconducting leads^ 
- a class of systems which plays a key role in the quest 
for Major ana fermions^^ - and ii) 2DEGs that form at 
interfaces between complex oxidea^l gpch as LaAlOa and 
SrTiOs, which display tunable SOCP^ and superconduc- 
tivit}El 

Motivated by this body of experimental and theoreti- 
cal literature, we investigate the collective spin dynam- 
ics of an archetypical 2DEG model Hamiltonian with 
Rashba SOC and s-wave pairing^, in the presence of re- 
pulsive electron-electron (e-e) interactions. In the ab- 
sence of superconductivity a Rashba 2DEG exhibits spin 
oscillations, which, at long wavelength and for weak re- 
pulsive interactions, have a frequency ^ 2aA:F, being 
the strength of SOC and the 2D Fermi wavenum- 
ber in the absence of SOC. These oscillations, however, 
are damped and quickly decay due to the emission of 
(double) electron-hole pairs, which, in the normal phase, 
are present at arbitrary low energies. In this Letter we 
demonstrate that in a Gor'kov-Rashba superconductor 
(GRSC), collective spin oscillations continue to exist in 
a wide range of parameters, and are undamped because 
they lie inside the superconducting gap where no other 
excitation exists. Fig. [l] shows schematically the nature 




FIG. 1: (color online) a) Response of a Cooper pair in the 
A = + chirality subband of a Gor'kov-Rashba superconductor 
subjected to an oscillating magnetic field in the y direction. 
The solid circle is the Fermi surface and the black dot is the 
origin of momentum space. The arrows labeled by "1", "2", 
and "0" describe the orientation of the spins under the action 
of a magnetic field that points up, down, or vanishes. Sponta- 
neous oscillations are sustained, in the absence of a magnetic 
field, by the internal exchange field, b) A supercurrent boosts 
the Fermi surface in the x direction (solid line) and creates 
a magnetic field in the y direction. As a result, spins begin 
to oscillate around the new equilibrium orientation, indicated 
by the thick red arrows. 



of the spin oscillations in a GRSC. At variance with the 
Cooper pairs of a standard s-wave semiconductor, the 
pairs of a GRSC are in a mixture of singlet and triplet 
states. It is this feature that enables the pairs to respond 
to an oscillating magnetic field applied, say, in the y di- 
rection. In the course of the oscillation the spins of a 
pair tilt in opposite directions, in a pair-breaking mo- 
tion that creates a net spin polarization along the y axis. 
The spin polarization produces an exchange field, which, 
if the electron-electron interaction is sufficiently strong, 
sustains oscillations of the appropriate frequency in the 
absence of an external field. The essential point is that 
these oscillations are undamped as long as their frequency 
falls below the quasiparticle gap: they will therefore dis- 
play an extraordinarily long lifetime^. 

In order to excite these long-lived spin modes one could 
in principle apply a short magnetic pulse, but there is 
also a purely-electrical method. Namely, a supercurrent 
pulse applied, say, in the x direction, will generate, via 
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the Edelstein effectP an effective magnetic field pulse in 
the y direction, and this should be sufficient to start the 
spin oscillations. This excitation mechanism is illustrated 
in Fig. ^p)' The Fourier spectrum of the supercurrent 
pulse must not contain frequencies of the order of (or 
larger than) twice the superconducting gap to avoid the 
creation of quasiparticle excitations. We suggest that the 
new collective spin mode can be used to realize "persis- 
tent spin oscillators" operating in the frequency range 
of 10 GHz — 1 THz (for superconductors with a critical 
temperature in the range 10~^ — 10 K). 

Model Hamiltonian and effective low- energy theory. — 
We consider the following model Hamiltonian: H = + 
+ y^e-e- Here V,o is the kinetic energy term given 
by % = V^J(^) hij{r) "ipjir), where {ft = I 

throughout this manuscript) 



hij{r) 
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- a [cTij X {-iVr)] -z- ji 5ij . (1) 



Here '0|(r) [^j(r)] creates (destroys) an electron with 
real-spin label i =t, | and band mass m, a measures the 
strength of Rashba SOC, cr = (<j^,(T^) is a 2D vector of 
2x2 Pauli matrices a^, /i is the chemical potential, and 
f is a unit vector normal to the 2D plane where electrons 
are confined to move (the x—y plane). Diagonalization of 
7^0 yields two bands, ^a(^) = h!^ / (2m)^\ak — ji^ A = ±1 
being the so-called "chirality" index. Rashba SOC forces 
spins to lie on the x — y plane and to be perpendicular 
to k at each point in momentum space [see Fig. [l^)]. 

The second term in the Hamiltonian Hp, is an s- 
wave pairing Hamiltonian which is responsible for super- 
conductivity: it physically corresponds to an attractive 
interaction of strength —g with ^ > 0, which is active 
only in a thin shell of momentum space around the Fermi 
surface. The microscopic mechanism responsible for the 
appearance of the pairing term is not important here. 
The problem defined by 1-Lq + Hp has been studied by 
Gor'kov and Rashba^ who calculated the in-plane and 
out-of-plane spin susceptibilities X||(±)(<? = ~^ 0)- 
Due to a mixture of spin-singlet and spin-triplet chan- 
nels stemming from SOC, the GRSC develops a finite 
and anisotropic spin response. 

In this Letter we study the spin response of a GRSC 
at finite frequency cj, taking into account also repulsive 
e-e interactions described by the last term in the Hamil- 
tonian H, 
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where V > and the spin-resolved density operator is de- 
fined by Pi{r) = 2pl{r)2pi{r). We are interested in study- 
ing the collective dynamics of the system described by 
V, assuming that it remains in a phase characterized by 
a hard (finite in any direction of space) gap, despite the 
presence of repulsive e-e interactions. These lead to an 
effective reduction of the parameter ^, in the spirit of the 
Anderson-Morel pseudopotential^^ . 



We now derive an effective low-energy action corre- 
sponding to the full Hamiltonian H in terms of spin 
degrees-of-freedom only. The first step is to decouple 
the two quartic terms. Hp and He-e, by means of a suit- 
able Hubbard-Stratonovich (HS) transformation (see e.g. 
Refs. I16|17|) . For the pairing term Hp we introduce the 
complex HS field Ao(r, r), which describes the supercon- 
ducting order parametei^^. We do the decoupling in the 
chiral basis: this allows us to work with Cooper pairs that 
are protected by time-reversal symmetr};^. Transforming 
back to the real-spin basis we get spin-triplet pairing in 
addition to the regular spin-singlet pairing^. 

It is useful to rewrite He-e as^. 



He 



V 



d^r 
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where p(r) = Pi{r) is the total-density operator, 
^a{'^) — j is usual spin-density op- 

erator, and C = {C1X2X3) is an arbitrary unit vector 
in 3D space. To decouple He-e by means of HS trans- 
formation we introduce four real HS fields^: ^(r, r) and 
M{r^ r), which are conjugate to density fiuctuations and 
spin fiuctuations, respectively. 

The notation is considerably simplified by defining a 
four-component spinor "^^{r^r) = [^ij ?/;| ip^] in real- 
spin space. The exact microscopic action corresponding 
to H after the HS transformation can now be expressed in 
a compact form as (the variables r,r will be suppressed 
from now on when needed for brevity) 



dr / d'r 
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where (3 = (A:bT)"\ So(r,r) = z0(r^ (g) 1^), and ^ is 
the Grassmann variable corresponding to the fermionic 
field Here Gq ^ is the Green's function of the problem 
defined by Ho + Hp^ and is a 4 x 4 matrix given by 



Go' = dtlr 



•r^ (g) + a {r X {-iV) • z} 
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The Pauli matrices act in the 2x2 Nambu-Gor'kov 
space and Icr (H-r) is the identity matrix in real-spin 
(Nambu-Gor'kov) space, T = {T^,T^,T^) = (r^ (g) 

(2) 

a^jtr g) cr^,r^ g) <j^) and A is a 2 x 2 matrix whose 
diagonal (off-diagonal) elements are related to the triplet 
(singlet) order parameter [see Eq. (A5)]. 

At low energies, fluctuations of the amplitude of the 
order parameter Ao(r,r) do not play any role while 
phase fluctuations give rise to the Bogoliubov- Anderson 
mod^. To this end, we write Ao(r,r) Ae^^(^'^\ with 
A real. The amplitude A is flxed by the saddle-point 
equation SS/ 5A = 0, which yields the BCS equatiorP 
[see Eq. ([A4|]. 



3 



The role of the phase field 0{r^r) can be made explicit 
in the action S by performing the following gauge trans- 
formation (fi{r^r) = ipi(r^T)e^^^'^'^^^'^ to new fermionic 
fields (fi{r^r). Writing the action S in terms of the new 
fermionic fields generates new self-energies in the round 
brackets in the second line of Eq. —Gq^ -\- J^q ^ 

Si + 112 + Ha with 
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Si(r,T) = 
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^2{r,T) = M T and E3(r,T) = ^[rx {Vr0)]-z . (7) 

The fermionic part of the action can be integrated out 
(since it corresponds to a Gaussian functional integral 
for the partition function) leaving us with the following 
effective action 
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FIG. 2: (color online) Panel a) - c) The in-plane dynami- 
cal spin susceptibility x\\{^j^) [i^i units of the 2D density- 
of-states m/{27v)] as a function of uj (in units of /i) for in- 
creasing values of A (in units of /x) and V = 0. The solid 
line represents ?R:e X||(0,cj), while the dashed line represents 
X||(0,<x;). Note that for finite A, X||(0,C(;) diverges at 
cj = cji and that X||(0,c<;) = for < cj < cji. Panel d) 
The quantity E-^{k)-\-E-{k) as a function of k (in units of ko). 
In this figure we have fixed a — 0.2/j/ko with ko = ^J^jnjl. 
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are the in-plane and out-of-plane dynamical spin suscep- 
Tr [in (— Gq + ^)] (S) tibilities of the GRSC described by Ho + '^p, respec- 



where the symbol "Tr" means a trace over all degrees of 
freedom (including space and imaginary time). 

To make further progress we need to expand the last 
term in tSeff in powers of S. We keep terms up to second 
order in the Fourier components of the fields (/>q,6>q and 
Mq. A remarkable simplification occurs in the q ^ 
limit where the action reduces to the sum of indepen- 
dent quadratic terms (see Appendix |B|. Density and su- 
percurrent oscillations on one hand and spin oscillations 
on the other hand decouple. As usual, the frequencies 
of collective modes are determined by the isolated poles 
of appropriate susceptibilities. For short range interac- 
tions, the density /current modes disperse linearly in q 
and their frequency vanishes at g = as expected for a 
regular Goldstone mode. The spin modes, on the other 
hand, have a finite frequency, which increases with in- 
creasing A [consistent with the fact that the resistance of 
Cooper pairs to the twisting motion described in Fig. [l^) 
increases with increasing A], but remains less than 2 A, 
ensuring long lifetime. 

Collective spin oscillations. — In the q ^ limit all 
the mixed response functions vanish (see Appendix [B| 
and the frequency of the collective spin mode c^n {oj±) at 
q = is given by the solution of the equation 



tively. These are obtained from the analytical continua- 
tion, iurn OJ -\- ^0+, of the corresponding expressions in 
imaginary frequency: 



Xa-ab{0,iUrn) 
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^Tr[r-Go(/c,i6, 
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where "Tr" implies a trace over spin and Nambu-Gor'kov 
indices and Vrn {^n) is a bosonic (fermionic) Matsubara 
frequency. After analytic continuation we find, at T = 0, 



X||(0,cj) 



and x±(0,cj) 
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: 2x||(0,a;) [£ = E^(k) ^ E_(k) and 
^yrv) = ^^^yrv ) ^ /S?]. Duc to the rclatlou between out-of- 
plane and in-plane spin response functions, we will dis- 
cuss only collective in-plane excitations. 



E\{k)^il{k) 
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X||(-L)(0,w) =0 



(9) 



with respect to uj. In passing, we note that Eq. ([9| can 
also be obtained diagrammatically from a vertex equation 
obtained by summing up ladder diagrams (see Appendix 
[C|. In Eq. ([9]), X|| = XaVi = Xa^a^ and x± = Xa^a^ 



We calculate x\\{^^^) numerically from Eq. (11) and 
plot its real and imaginary parts in Fig. [2] In tne limit 
A = {i.e. absence of superconductivity) - see panel 
a) - the imaginary part is non-zero only in the interval 
of frequencies between 2aA:F,+ and 2a/cF,- [^f,± being 
the minority (majority) Fermi wave vectors for the two 
Rashba bands ^a(^)] and the real-part exhibits (logarith- 
mic) singularities at these boundaries (see Appendix [D|). 
When this result is inserted in Eq. ([9|, one finds a col- 
lective spin mode, which is undamped within this ap- 
proximation. In a more refined theory (beyond Gaussian 
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operating in the frequency range of 10 GHz - 1 THz. 
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FIG. 3: (Color online) Panel a) The cj ^ limit of X||(0,^) 
as a function of a/co/A. The solid line is the result obtained 
from Eq. (11) while the filled circles are the result of Ref. [G] 
Panel b) A 2D color plot of the frequency c<;|| of the in-plane 
collective spin mode (in units of A) as a function of the in- 
verse of the strength of electron-electron repulsions (Vb/V^, 
with Vo = 47r/m) and SOC (aA:o/A). In this plot A = 0.1 \i. 
The top contour line is for l>j\\ — 2 A while the bottom con- 
tour line defines the boundary of the region in which =0. 
The collective spin mode is undamped when it lies within the 
superconducting gap (0 < L)J\\ < 2A), i.e. when falls in the 
region enclosed by the two contour lines. Panels c) and d) 
represents ID cuts of the plot in panel b). 



fluctuations), however, low-energy double electron-hole 
excitations damp this mode. We now show that, at odds 
with the normal phase, in the superconducting state the 
mode lies (for a wide range of parameters) within the su- 
perconducting gap and thus cannot be damped by these 
excitations. 

In panels b) - c) we plot x\\{^^^) fo^ finite A. In 
the superconducting state 5Re x\\{^^^) exhibits a diver- 
gence at ui = min/c[£^+(/c) + E_{k)]. In panel d) we 
plot Ej^{k) + E-{k) as a function of k. In the region 
< a; < cji, x\\{^^^) is identically zero and, since 
5Re x\\{^^^) diverges for ^ c^i, there is always an in- 
plane collective spin mode with frequency uj\\ ^ Ui for 
weak repulsive interactions V. Our results for the fre- 
quency of the in-plane collective mode uj\\ as a function of 
V and a (for a fixed value of A) are summarized in Fig.js] 
Note that there is a wide range of parameters such that 
lies within the superconducting gap, < < 2A. 
We also have checked that, as expected, uj\\ increases with 
A. 

In summary, we have shown that quasi-two- 
dimensional superconductors with tunable spin-orbit 
coupling exhibit undamped collective spin oscillations 
that can be excited by the application of a magnetic field 
or a supercurrent. The concerted action of spin-orbit 
coupling and electron-electron interaction is essential to 
the establishment of these collective oscillations. Since 
the frequency uj\\ of these oscillations is of the order of 
the superconducting gap A we expect that our findings 
might enable the realization of long-lived spin oscillators 
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Appendix A: The Green's function of a 
Gor'kov-Rashba superconductor 

Let us first consider the Green's function Gc of a 
Gor'kov-Rashba superconductor in the so-called "chi- 
ral" basis. We remind the reader that we define 
"Gor'kov-Rashba superconductor" the system defined by 
the Hamiltonian Ho-\-Hp^ where the Rashba Ho and pair- 
ing Hp Hamiltonians have been defined in the main text. 

We start by defining the four-spinor ^J(/c) = 
[ip\.{k) 2p-^{—k) 2p^_{k) 2p-{—k)] in momentum space, 
where ipx{k) and ip\{k) are field operators correspond- 
ing to the eigenstates of the Rashba Hamiltonian Ho 
introduced in Eq. (1) of the main text. In this "chi- 
ral" basis the Matsubara Green's function correspond- 
ing to the Gor'kov-Rashba Hamiltonian Ho + Hp [i.e. 
Gc{k,r) = -{Tr{'^^{k,r)^{k,0)})] is a 4 x 4 block- 
diagonal matrix (in Fourier transform with respect to 
imaginary time r): 



G^{k,ien) 

G_ 





{k.icn) 



where Gx{k,ien) are the following 2x2 matrices: 
-icntr - ^x{k)r^ + AAr • k/k 



Gxik.icn) 



■El{k) 



(Al) 



(A2) 



Here El(k) =^l{k)^A^ and Cn is a fermionic Matsubara 
frequency. The gap A of the Gor'kov-Rashba supercon- 
ducting state is given by 



(A3) 



where ^ > is the pairing coupling constant (see 
main text) and (pk is the angle between k and the 
X axis. We emphasize that, following Ref. |6j we 
have assumed that pairing occurs only between time- 
reversed partners within each Rashba spin-orbit-split 
band: ^)'0a(^)) = if A = —A. The gap A is fixed 

by the saddle-point equation SS/SA = 0, [5 is given in 
Eq. (4) of the main text] which yields the mean-field BCS 
equation^ 
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k,X 



tanh [(3Ex{k)/2] 
'2Ex{k) 



(A4) 
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where A is the system's area. 

The Green's function Go introduced in Eq. (5) of the 



main text is exphcitly given by the following 4x4 matrix: 



/dr-Vl/{2m) a{d^-idy) An A12 \ 

-a{d^ + idy) dr-Vl/{2m) -A12 A22 

An A12 dr + Vl/{2m) a{d^ + idy) 

V -A12 A22 -a{d,-idy) dr + Vll{2m)J 



(A5) 



where Aii(r,r) = (?/^^(r, r)?/^^(r, r)) and A22(^) = 
{%!) ^{r ^ t)iI) ^{r ^ t)) are related to the triplet order pa- 
rameter (which, of course, arises only because of 
the presence of spin-orbit coupling) and Ai2(r,r) = 
('0^(r, r)'0>^(r, r)) is the singlet order parameter. 

The 4x4 Green's function Go in the real-spin basis 
can be related to the Green's function in the chiral 
basis: we find (in Fourier transform with respect to space 
and imaginary time) 



Gl^{k,ien) = F^{k,ien)(j'^-F^(k,ien) (k x - z 
and 

Go^(/c, icn) = Gs{k, ien)la - G^{k, icn) (jsx a 
Here k = k/\k\ and we have defined 



(A6) 
(A7) 

T 

(A8) 



(A9) 



2G,/a(fc, ien) = (k, ie„) ± G]} (k, ie„) (AlO) 
-ien - , -ien -^-(fc) 



2G,/^{k,ie„) = Gf{k,i€„) ±G^_^{k,i€„) (All) 

^ -ten , -Kn+£.-{k) 

" e2+e+(fc)2 + A2 £2 +e_(fc)2 + A2 ' 



and 



2F,/,{k,ien) = e'^^ [G^_^{k,ien) ^ G^^{k,ien)\ (A12) 

1 



4 + e-(^P + A2 



In Eqs. (A6|-( |A12| Gq^ (G±^) are the elements of the 
matrix Gq (G±)'. 
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Appendix B: The effective action, Gaussian 
fluctuations, and spin coUective excitations in the 
long- wavelength limit 

To include the phase fluctuations of the order parame- 
ter in our study (ignoring the fluctuations in the modulus 
of the order parameter) , we perform a gauge transforma- 
tion to new fermionic fields 



and 



(^l(r,r) = ^|(r,r)e-^^("'")/2 , 
(^,(r,r) = Vi,(r,r)e^^^"'")/2 . 



(Bl) 



(B2) 



Using Eqs. (|B1|)-(B2) in the definition of the action S 



given in Eq. (4) of the main text we find: 



S 
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:(-Goi+Ei + S2+S3) 



(B3) 



where (?(r, r) = [(^| (^| (f)^ (p^ is a four component spinor, 
and Si, 112 and are given by Eqs. (6) - (8) of the 
main text. Now the fermionic fields (^) can be inte- 
grated out of the partition function, which is of the form 
Z = Zq j D[^] exp {3B<P) /2, by performing a Gaussian 
integral over the Grassman variables. We use the follow- 
ing relations for an arbitrary matrix 



j D[^]exp Q^^^^ = \/Det[5] 



exp 



■ln(Det[5]) 



= exp ( -Trln[5] ) . (B4) 



Using Eq. (B4) in Eq. (B3) immediately gives the effec- 
tive action 5eff reported in Eq. ([8| of the main text. 

To expand 5eff up to second order in S, we use the 
identity 



Tr [In (-Go^ + S)] = Tr [in {-G^')] 



Tr 



OO -J 



(B5) 
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Using Eq. (8) in the main text and the identity above we 
find, up to second order in S, 



->eff 



drdr 
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ixr [In (-Go^)] + ^TV[GoE] + ilV[GoEGoE] 



(B6) 



where the trace "Tr" is taken over ah degrees of freedom 
including space and imaginary time. Note that in the 
last term of the previous equation we have to keep only 
terms up to second order in the fluctuating fields 6>, 0, 
and M (Gaussian fluctuations). 
We can rewrite <Seff as 



5eff ^ j drdr^ - Tr [in (-^o ')] + 



(2) 

fluct ' 



(B7) 



(2) 

where ^S^^^t second-order contribution to 5eff 

due to the fluctuating fields. Using Fourier transforms 
with respect to space and imaginary time and defining 
Ii{q,ii^m) = [0{q,ivm) (l){q,ii^m) M{q,iVrn)] we find 

4uct = ^n(gr,i^^) N{q,iu^) IV^{-q,-iym) , (B8) 

q,m 

where A/" is a 5 x 5 matrix whose elements are given by 
various dynamical response functions and Vrn is a bosonic 
Matsubara frequency. Different elements of the matrix 
N{q,iUm) are given by 



Pmf 

m 



-Xrjj{q^i^m) 



N22{q,i^m) = ^ + ^Xppiq^i^m) 



(BIO) 



Ni2{q,iiym) -^^mXppiq^^^m) 



(Bll) 



-^^^QiX^pi-q^-i^m) + ^<7X0p(-<?, -^^m) , 



N2l{q,iiym) = -^^mXppiq^^^m) 
Nab{q,ij^m) = ^^ab - ]pCa--2 ^b-2{q, IVrn) , (B13) 



with a, 6 G [3, 4, 5], 



Nia{q,i^m) = -^^mXa-p(9'^^m) 



^QiXa-j-iq^^^m) - jqXa-c/>{q^ij^m) , (B14) 



Nal{q, i^m) = ^^mXa— 2p(~^' -^^m) (B15) 
--lY^iXcr — ^j-{-q,-^^m) - ^qXcr—^c^{-q,-^^m) , 



N2a{q,i^m) = -iXa— ^ (B16) 

and, finally, 

Na2{q, i^m) = -^Xa«-2p(-g, -^i^m)/2 . (B17) 

The quantity pmf which appears in Eq. ( |B9[ ) is given by 
= ^ Tr[^o(/^, en)r^ la] , 



(BIS) 



and physically corresponds to the total electron density 
(superfluid and normal component^. 



Th e response functions x that appear in Eqs. (B9) 
(B17) are dynamical susceptibilities of the Gor'kov 



Rashba superconducting state in the absence of electron- 
electron interactions: 



Xpp{q.^^m) = -^^Tr [Go(/c+,e+) ^ 1, Go(/c",e-) 1, ] 



k,n 



(B19) 



Xpp{q,ii^m) = -Ajl^^ Tr [Go(fc+,e+) 1,(^1^ Go{k-,e-) 1, ® 1^ ] , 

^ k,n 



(B20) 
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1 X k- 

X...4g,^^^m) = -^E^r [Go(fc+,e+)P Go(fe-,e-)r] , 
Xav(9' ^'^'") = E [Go{k+, e+) t3 $5 1, Go(fc-, e-)r] , 
X<,.,-.(q,ir^™) = - ^ Tr [Go(fc+,e+) 1, 65 1, Go(fe-,e;;)r] , 



Mq,w^) = -^E^r [Go(fc+,e+)r^ Go(fe-,e;;)P] 



{q,iv^) = E'l^ [Go(fe+,e+) t3 ® 1, Go(fc-,e-) r<^] 



k,n 



and, finally, 



^(g,it^m) = -^E^ Tr[Go(fe+,e+) 1,01, Go(fe-,e-)r^] 

k " 



(B21) 
(B22) 
(B23) 
(B24) 
(B25) 
(B26) 
(B27) 

(B28) 



r 



In Eqs. (B19)-(B28) fc^ = fc ± q/2, = e„ ± i^mf^, the ductor in the absence of spin-orbit-coupling (a = 0), the 
trace "Tr" is to be taken over the Nambu-Gor'kov and quantities Nii,Ni2,N2i,N22 reduce to the elements of 



spin indices, and the 4x4 matrix is 



/ e'^^ \ 

_e-^0«j 

e-^^^ 

\ -e*^^ J 



(B29) 



the matrix B{q) defined in Appendix A of Ref.lTTland de- 
termine the frequency of the Bogoliubov- Anderson mode 
once the limit q ^ is taken. 

To find the energy-momentum dispersion of the col- 
lective spin excitations, we obtain an effective action in 
terms of the spin degrees-of-freedom M only by per- 
forming a. Gaussian integral over the fields (f) and 



pressed in a compact form, by defining lV~{q^iUrr 
[Mi{q,iUm) M2{q,ivm) Ms{q,iUm)]^ as 



The matrices with i = 1...3 have been introduced in Eq. (B8). The "spin-only" action S)J can be ex- 
in the main text. Notice that the 4x4 matrix (8) l^r 
corresponds to the density operator, the 4x4 matrices 
correspond to the spin operators 5^, and the 4x4 matrix 
(It- (g) lcr)ki/m corresponds to the current operator ji. 

In deriving Eq. ( |B8| ) we have used the following rela- 
tions: XAB{q,i^m) = XBA{-q, -i^m) whcrc Aoi B de- 
note density, spin or current operators and X(f)A{Q^ i^m) = 
—XA(f){—Q^ —i^m)' Note that for an ordinary supercon- where Q is a 3 x 3 matrix whose elements are given by 

I 



q,'m 



(B30) 



Qcd{q, iVm) = iVc+2 d+2 + [Ni d+2{N22Nc+2 1 - N2lN^+2 2) + N2 d+2{NnNe+2 2 - N12NC+2 1)] , (B31) 



with c,d E [1,2,3] and where A{q,iVm) = N11N22 — -^2i-^i2- Note that the first term in Qcd originates from 
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the correction to the bare dynamical response function 
due to electron-electron interactions, while the second 
term originates from the coupling of spin fluctuations to 
phase fluctuations. The collective spin modes of the sys- 
tem can be found by solving 



(B32) 



As mentioned in the main text, in this work we are 
interested in finding the frequency of the collective spin 
mod es in the long- wavelength q ^ limit. In this limit 
Eq. ( |B8[ ) simplifies considerably and the effective action 
decouples into separate terms corresponding to super- 
current /density oscillations and spin oscillations, respec- 
tively. More specifically in the q ^ limit we have. 



(B33) 



Moreover, 



Xa^a^ (gr, iUrn) = Xa^piQi ^^m) = (B34) 

for every finite q. In the limit q ^ Eqs. (B33)-( B34| ) 



give 



and 



= 7V24 = A^25 = Ns2 = N^2 = A^52 = , (B35) 



(B36) 



A^34 = iV35 = iV43 = A^45 = N^?> = A^35 = . 

Using Eq. ( |B35D in Eq. \mi\ we obtain 

Qcd(0, ivm) = Nc+2 (i+2(0, ivm) • (B37) 

In other words, phase/density fluctuations do not couple 
to s pin flu ctuatio ns in the long- wavelength limit. Using 
Eq. ( |B36[ ) in Eq. ( |B31| ) we obtain Q^d = for c 7^ c^, i.e., 
all the off-diagonal components of the matrix Q are zero 
in the long wavelength limit. Eq. (B32) thus reduces to 



det[Q(0, ivm^oo^ iO+)] = 

Xa2a2(0,Cj) 



2 
V 



2 
V 



2 
V 



Xa3a3(0,Cj) 



= , (B38) 



with X(Tiai(0,Cj) = Xa2a2(0,Cj). 

Note that we ca n also obtain Eq. (B38) directly from 

(2) J I 

<S^^^^^^ Equation (B35) implies that the matrix N in 
Eq. (B8) has a block diagonal form comprising an up- 
per 2x2 block corresponding io — (j) fields and a lower 
3x3 block corresponding to the M fields. Thus the 
q = {) component of the action S^^^^ can be expressed 
as a product of a "phase only" action and a "spin only" 
action, i.e. 



c(2) 
^fluct 



- C(2) 



g=0 



g=0 



X S 



q=0 



(B39) 



Moreover, Eq. (B36) implies that the lower 3x3 block of 



the N matrix in Eq. (B8) is diagonal in the long wave- 
length limit. We thus obtain 



S 



q=0 



2 
V 



XM^(0, iUrr 



(B40) 



which gives the same condition for the existence of col- 
lective spin modes as Eq. (B38). 



Appendix C: The ladder sum and the vertex 
equation 

In this Section we show that the equation 

2 2 

y - Xai^i (0, uj) = — - XaV2 (0, uj) 



V 



Xa3a3(0,Cj) = 



(CI) 



for the collective spin excitations that we found in the 



previous section, Eq. (B38), can also be obtained dia- 
grammatically. 

In the (conserving) ladder approximatioiP^^ the dy- 
namical spin response function Xa'^a'^ the presence of 
electron-electron interactions is given by 

Xa^a^iq.ii^m) = - ^ ^ Tr ^E^Gq (A^+ , ^6+) A(gr, zi^^) 



Go(fe ,ie„)j . 



(C2) 



The vertex function K{q,iVrn) is a "dressed" version of 
the bare vertex = 1^- (g) and it accounts for the 
interplay between electron-electron interactions and the 
external electromagnetic field. Similar equations hold for 



Xo 



and Xc 



The vertex A is a 4 x 4 matrix and satisfies the following 
equation (see Fig. |4|: 

K{q,iVm) = r2-FT3®l,Jl^Go(fc+,ie+) 



X A(gr,zi/^)Go(fe",ie-)| r^0l^ . (C3) 

In the g ^ limi t, af ter some lengthy but straightfor- 
ward algebra, Eq. (C3) yields 



K{{),iVrn) 



V 

1 - ^XaV2(0,ii^m) 



(C4) 



In the ladder approximation, the interacting in-plane 
spin-susceptibility in the long-wavelength limit is thus 
given by 



Xa^a^iO.iUm) 



Xa^a^{0,iUrr 



V 

1 - ^Xa2a2(0,ii^m) 



(C5) 
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k + q/2 



k + q/2 




k-q/2 



k-q/2 



FIG. 4: Panel a) The bare bubble Xa'^a'^ cons tructe d with the 
Gor'kov-Rashba Green's function Go [see Eq. (B22 )]. Panel b) 
The dynamical response function Xa'^a'^ with the inclusion of 
vertex corrections due to electron-electron interactions. Panel 
c) The vertex equation in the ladder approximation. In this 
figure, k ± q/2 denotes the 4 vector {cn ± i^m/'^, k =b q/2) and 
= Ur cr^ denotes the spin operator §2. 



Since collective modes are isolated poles in the dynami- 
cal response function ^^2^2 (0, zi^rn OJ -\- iO~^) (located 
infinitesimally below the real- frequency axis), Eq. ( |C5| ) 
reproduces the condition given above in Eq. (|C1[). 



Appendix D: Normal Rashba Gas 

In this Section we report explicit expressions for 
the real and imaginary parts of the in-plane dynam- 



ical spin susceptibility Xcr2cr2 (0, a;)|^^Q of a normal 
(non-superconducting) Rashba 2DEG in the absence of 
electron-electron interactions [see panel a) of Fig. 1 in the 
main text]. 

At zero temperature and in the absence of supercon- 
ductivity we find: 



A=0 



1 r^^'- 



kdk 



1 



CO - 2ak + iO+ 



uj + 2ak + zO+ 

(Dl) 



where /cf,± = (2m/i + m^a^)^/^ ^ ma is the Fermi wave- 
vector of the minority (majority) Rashba band. Perform- 
ing the integration in Eq. (Dl), we find that the real and 



imaginary parts of Xa^a^{^y^)\A=o given by 



m 



A=0 



1 



X In 



27T I ' 8ma^ 



and 



UJ — 2a/cF,+ UJ + 2aA:F,- 

(D2) 

6(0; - 2aA:F,+) 

uj) . (D3) 



^=0 16^2 
X 6(2aA:F 



Eqs. JD2| and (|D3]) agree with Eqs. (7) and (10) in 
Ref. [2Q,Tafter setting to zero the Dresselhaus spin-orbit 
coupling constant in their results). 
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